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Abstract 

We apply an integral inequality to obtain a rigorous apriori estimate of the accuracy 
of the partial sum to the power series solution of the celebrated Riccati-Bernoulli 
differential equation. 

1 Introduction 

The deservedly celebrated Riccati-Bernoulli differential equation: 



/ 2 2 

y = X +y 



XI) 



has the general solution ([I]): 



y[x) 



J3 (ix^) - cJ^ {\x-) 
'cJ.{\x-)+J^{\x-) 



(1.2) 



where c is an arbitrary constant and where 



JJx) :-- 



X 



2"r(n + 1) 



X 



X 



22.l!.(n+l) 24.2!-(n + l)(n + 2) 



;i.3) 



is the Bessel function of the first kind of order n, where n is any real number, and 
r(n + 1) is the famous gamma function. The very interesting history of (11. ip is detailed 
in Watson's standard treatise [3J. 

The equation (11.11) is an example of a simple differential equation whose solutions form 
a family of transcendental functions which are essentially distinct from the elementary tran- 
scendents. 
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Unfortunately, the general solution (11.21) does not easily lend itself to a rigorous error 
analysis of its accuracy in a particular interval of the variable. 

We will show how a simple application of an integral inequality allows one to estimate 
the accuracy of the partial sum of the Taylor series expansion of the solution within the 
latter's interval of convergence. 



2 Cauchy's theorem 



(In this section we follow |2], Chapter IV, section 5). The general Cauchy Problem is to 
solve the ordinary differential equation (ODE) initial-value problem: 



y' = fix,y), y{xo):=yQ 



Let f{x,y) be expanded in the series 



fix, y) = ^ Ajix - xoYiy - yoY 



convergent for 



\y - yo\ < R2 



(i?l > 0, i?2 > 0) 



(2.1) 



(2.2) 



(2.3) 



Then, according to Cauchy's theorem in the theory of differential equations, the problem 
(12. ip has a solution y{x) represented by the series 



(2.4) 



convergent in some neighborhood of the point Xq. 

Cauchy's theorem allows one also to indicate the neighborhood of the point xq, in which 
the series (12.41) converges. Namely, let M be a constant such that 




\f{x,y)\^M 



with 



\x — XqI ^ Ti < Ri 



\y - yo\ ^ rs < i?2 



(2.5) 
(2.6) 



where Ri and R2 are numbers, defining the region (12.31) of the convergence of the series 
and ri and r2 are some positive numbers. Then the series (12.41) converges for 



\x — Xo\ < r 



where 



ri 1 



(2.7) 



(2.8) 



It should be noted that the true interval of convergence is usually much larger than (12.81) . 
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3 An integral inequality 

We use the notation of (12. ip : 

Theorem 1. Let I be the interval, xq ^ x ^ xi, and suppose that for all x E I, f{x, y) > 0, 
and that the differential inequality 



y' ^ f{xi,y{x)} 
also holds there. Then, the integral inequality 



Xo 



dy 



fixuVit)] dt 



dt ^ x — Xo 



{3.1] 



(3.2) 



holds for all X E I. 

Proof. This is a simple consequence of the elementary calculus sufficient condition that a 
function be decreasing in an interval. 
Define 

Taking the derivative and using the fundamental theorem of calculus we obtain 

1 dy' 



9 [X) 



But, the inequality (13. ip shows that 



for all X E I. 



f[xi,y{x)] dx 



g\x) ^ 



(3.4) 

(3.5) 
□ 



4 The Riccati-Bernoulli Initial Value Problem 

(In this section we follow [2], Chapter IV, section 5, but with some important refinements.) 

Problem: It is required to find the first 11 terms of the power series expansion of the 
solution of the initial value problem 



, 2 , y 
y = X + — 



y{0) := -1 



(4.1) 



and an interval of convergence for it. If we specialize the general solution (11.21) to this initial 
value problem we find the following formula for the exact solution: 



y{x) 



^ f r(!)J|(|x^)-v^r(i)j^(ix^) 
16 ■ 1 v/2r (i) J. {\x^) + r (I) {\x^) 



(4.2) 
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Unfortunately, the computation of the power series solution on the basis of the quotient 
fl4.2p . although theoretically possible, is computationally formidable. 

It is easier to use the equation fl4.ip to compute the derivatives of y{x) at x = directly: 

y'= x' + i =1 (4-3 

y"= '^x+\yy' (4.4 

67 

y"'= 2 + \y'' + \yy" =^ (4.5 

2/^'^= ly'y" + \yy"' = (4.6 

207 

^(5) = yi2 + 2y'y"' + ly^/W _ ^ 7 



,(6) 



128 

5^Y" + §2/ V^) + =-—r (4.^ 



64 

y(^) = + f + 3y'y(^) + iyy(^) = ^ (4.9 

= f + f + ly'yi^) + = (4.10 

^(9)^ f(y(4))2 + 28/V^) + WV^)+4yV^) + |y?/W =^|^ (4.11 

^(10) = 63?/Wy(^) + 42y'>(6) + 18y"y(^) + + = -^^^ (4.12 
Therefore, by the formula fl2.4p . the first 11 terms of the series solution are: 



^ 1 1 2 67 3 35 4 69 . 77 6 5317 . 

y[x) = — 1 H — X X H X X H x x H x 

' 4 16 192 768 5120 15360 1032192 

2655 „ 80641 „ 77851 

-x^ H —x^ H 



(4.13) 



1835008 165150720 62914560 



To find an interval of convergence of this series (14.131) we use the Cauchy theorem. It x 
and y satisfy the inequalities 

|a;| ^ 0.5 |2/ + 1|^1 

then we may conclude that 

\f{x,y)\ ^|0.25[(y + l)-l]2 + x2| 
^0.25(1?/ + 1| + lf + \x\^ 
^1.25 

Therefore, in the formula (12. 8p . we may take 

ri:=0.5 r2 := 1 M := 1.25 
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and the value of r we obtain is: 

r = 0.5 (1 - e-°-^) = 0.2753355 ■ ■ • (4.14) 
Therefore, the power series solution (14.131) most certainly converges for \x\ ^ 0.27. 

5 The Accuracy of a Partial Sum from the Integral Inequality 

We will consider the following concrete problem although the principles are of general appli- 
cability. 

Problem: It is required to determine the accuracy of the partial sum of degree 9 of the 
power series solution ( 14.13P in the interval ^ x ^ 0.2. 

Since the series (14.131) is the Maclaurin expansion of y{x), we must estimate the re- 
mainder term, Rg{x), which we write in the LAGRANGE form: 



Rq(x) := ; X 

^ 10! 



(5.1) 



where < ©g < 0.2. We have to estimate y^^'^\x), the formula for which is given in (I4.12p . 
for all values of x in the interval ^ x ^ 0.2. The formulas (14. 3 p through (14. lip show us, 
finally, that we must estimate y{x), itself in ^ x ^ 0.2. 

The estimate of y{x) via our integral inequality (13. 2p constitutes the novelty in this paper. 

Maintaining the notation of (13.21) we see that 

Xq ■= xi := 0.2, (5.2) 

that the right hand side 

fix. n) := x'^ -\- ^ 

on /, and the differential inequality (13. ip becomes: 



fix,y):=x' + ^>Q (5.3) 



dx 4 
Therefore, the integral inequality (13.20 becomes 

1 dy 



^<0.04+^ (5.4) 



n I y(ty dt 

° ' 0.04+ ' 



dt^x (5.5) 
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But, 



' ^^Oarctan^^ (5.6) 



, y(tf dt { dt { 0.4 
0.04 + ' 

4 



5 



Therefore, 

y(x) —1 

10 arctan 10 arctan — ^ x 

0.4 0.4 

or 

y{t) X —1 

arctan < h arctan — - 

0.4 10 0.4 

and taking the tangent of both sides and reducing we obtain the estimate 




(5.7) 
(5.8) 

(5.9) 



which holds for all x G /, 

The function on the right-hand side of (15. 9p is monotonically increasing in /, and 



tan (f ) - 1 



1 + 



0.2 



-0.9447608 ■ ■ ■ < -0.94 



and we have proved that the following inequality is true for all x E I: 



-I ^ y{x) ^ -0.94. 



(5.10) 



{Note: the true value of 2/(0.2) is 



2/(0.2) = -0.9497771- •• , 
so the estimate (15. 9p . with an error of —0.00501 ■ ■ ■ , or about 0.53%, is quite good!). 



Now we must estimate y'{x) 



y 



(10), 



x) using the formulas (14.31) through (14.121) . 



Theorem 2. The following inequalities are valid for all x E I. 



-1 


^ y{x) 




-0.94, 


0.22 


^ y'ix) 




0.29, 


-0.15 


^ y"ix) 




0.3, 


1.87 


<: y"\x) 




2.12 


-1.13 


^ y^^Kx) 




-0.74, 


1.17 


^ y^^\x) 




1.93, 


-3.38 


^ y^^\x) 




2.23, 


14.59 


^y^'\x) 




27.12, 


-61.96 


^y^'\x) 




-22.73 


92.03 


^y^^Kx) 




146.76, 



and finally 



-665.9 < y^^^\x) < 281. 



(5.11) 
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Proof. All of our estimates come from worst case values applied to each of the summands in 
the formulas. 
By (m 



2 , y 



Therefore, using flS.lOp and (15.1 21) . we conclude that for all a; G /, 

y' < 0.2^ + 



(5.12) 



while 

y'>o' 

Therefore, we obtain the bounds 



-.94)^ 



.2209 > .22. 



.22 < y'{x) < .29. 



By (m 



If 1 / 

y =2x+ -yy . 



Therefore, using (I5.13P and (15.141) . we conclude that for all x G /, 



while 



y" < 2(0.02) + ^(-0.94)(.22) = 0.2966 < .3, 



y" > 2(0) + ^(-1)(.29) = -0.145 > -0.15. 



Therefore, we obtain the bounds 



By (m 



-0.15 < y"{x) < 0.3. 



/// o I '2 I // 

y =2+2?/ + • 



Therefore, using (I5.13P , (15.150 and (15.160 . we conclude that for all x E I, 



while 



y'" < 2 + ^(0.29)2 + i(-0.15)(-l) = 2.11705 < 2.12, 



1 1 

y'" > 2 + -(0.22)2 _^ _(o.3)(-l) = 1.8742 > 1.87. 

2 2 



Therefore, we obtain the bounds 



(5.13) 



(5.14) 



(5.15) 



(5.16) 



1.87 < y"'{x) < 2.12. 



(5.17) 



By gS]) 



- ^yy +^yy ■ 



(5.18) 



Therefore, using fl5.13p , (15.151) . (I5.17p . and (15.181) . we conclude that for all x E I, 

< ^(0.29)(0.3) + ^(-0.94)(1.87) = -0.7484 < -.074, 



while 



y(^) > ^(0.29)(-0.15) + ^(-1)(2.12) = -1.12525 > -1.13. 



Therefore, we obtain the bounds 



-1.13 < y(^)(a;) < -0.74. 



By (SZ!) 



y^''' = \y"'^^yy"^\yy^'y 



(5.19) 



(5.20) 



Therefore, using (157[3|) , (15351) . (ISTTj) . (15191) and ([E2QD, we conclude that for all x G /, 
?/(^) < ^(0.3)2 + 2(0.29)(2.12) + ^(-1)(-1.13) = 1.9296 < 1.93, 



while 



^(5) > ^(0)2 + 2(0.22)(1.87) + ^(-0.94)(-0.74) = 1.1706 > 1.17. 



Therefore, we obtain the bounds 



1.17 < y^^\x) < 1.93. 



(5.21) 



By 



(5.22) 



Therefore, using (l5l3D , (15151) . (l5Tfl) . (l519l) . (lOTD and (15:221) . we conclude that for all 
a; e /, 

< 5(.3)(2.12) + ^(.22)(-0.74) + i(-0.94)(1.17)) = 2.2231 < 2.23, 

while 



y(^) > 5(-0.15)(2.12) + -(.29)(-1.13) + -(-1)(1.93)) = -3.37425 > -3.38. 



Therefore, we obtain the bounds 



3.38 < y^^\x) < 2.23. 



(5.23) 



By flQ]) 



(5.24) 



Therefore, using (15^^ , flET^ . flCTTD . flEm . fICT]) . flHT^ and flOD . we conclude that 
for all X E I, 

?/(^) < 5(2.12)2 ^ ^(_o.i5)(_i.i3) + 3(.29)(1.93) + i(-l)(-3.38), 



which 
while 

which 



27.11235 < 27.12, 



y^'^ > 5(1.87)2 + ^(0.3)(-1.13) + 3(.22)(1.17) + ^(-1)(2.23), 



14.5992 > 14.59. 



Therefore, we obtain the bounds 



14.59 < y^'^\x) < 27.12. 



By glOD 



(5.25) 



(5.26) 



Therefore, using flSTTSj) , fl57[5|) . flSTTj) . flSTTgl) . flOB . dESS]), fl5:25|) and flOGl) . we conclude 
that for all x G /, 

35 21 7 1 

yis) < _(i.87)(_o.74) + _(o.3)(1.93) + -(.29)(2.23) + -(-0.94)(14.59), 



which 
while 

which 



-22.73085 < -22.73, 



35 91 7 1 

> :^(2.12)(-1.13) + -(-0.15)(1.93) + ^(.29)(-3.38) + -(-1)(27.12), 



-61.95345 > -61.96. 



Therefore, we obtain the bounds 



-61.96 < y^^\x) < -22.73. 



(5.27) 
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By gUD 



(5.2^ 



Therefore, using (l57[3D , (15151) . (ISTTD . (15191) . (ICT]) . (15123]) . (lOSj) . (151271) and (15:281) . we 
conclude that for all x E I, 



yi9) ^ + 28(2.12)(1.93) + 14(0.3)(2.23) + 4(.29)(27.12) + -(-1)(-61.96), 



35 



which 
while 

,{9) 



146.75575 < 146.76, 



35 



r ^ > ir(-0-74)^ + 28(1.87)(1.17) + 14(-0.15)(2.23) + 4(.22)(14.59) + -(-0.94)(-27.73), 



which 

Therefore, we obtain the bounds 
By gl2D 



92.0335 > 92.03. 



92.03 < y^^\x) < 146.76. 



(4)„{5) 



+ 42/'2/(6) + 18/2/(^) + -y'y^^^ + -yy^'\ 



(5.29) 



(5.30) 



Therefore, using ^J3i) , ^M, (M), ^M, ^M, ^M, ^M, and (00]) . 

we conclude that for all x E I, 

?/(^°) < 63(-0.74)(1.17)+42(2.12)(2.23)+18(0.3)(27.12) + ^(0.22)(-22.73) + ^(-0.94)(92.03), 



which 



while 



280.9922 < 281, 



y^^°^ > 63(-1.13)(1.93)+42(2.12)(-3.38)+18(-0.15)(27.12)+^(0.29)(-61.96)+^(-l)(146.76) 
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which 



Therefore, we obtain the bounds 



-665.8137 > -665.9. 



-665.9 < y^^^\x) < 281. 



The inequality (I5.3ip was the goal is this long and detailed computation! 



(5.31) 



□ 
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Now we can state the accuracy of the partial sum: 
Theorem 3. For all x in ^ x ^ 0.2, the partial sum of degree 9: 

, 1 1 2 67 o 35 4 69 . 77 e 5317 . 

y[x) «i — 1 H — X X H X X H x x H x 

' 4 16 192 768 5120 15360 1032192 /c- oo^ 

2655 8 80641 g ^^'^^ 
~ 1835008^ ^ 165150720^ 

approximates the true value, y{x), of the solution series, (14.131) . with an error that does not 
exceed 2 units in the eleventh decimal place. 



10! 



x'' 



^ ^(0.2)10 = 1.878 ■ • ■ 10-" < 2 ■ 10-" 



□ 



Proof. The estimate 

\R9ix)\ = 

completes the proof. 
6 Conclusions 

Our integral inequality (13.21) can be applied to wide classes of differential equations. 

For example, our method allows us to prove that in the interval ^ x ^ 0.4 the polyno- 
mial 

Vix) := IH \ H x^ H x^ H x^ 

' 4 16 192 96 200 

approximates the true solution, of the initial value problem 

4y' = a; + y^ 1/(0) :=1 

with an error that does not exceed 2 units in the fifth decimal place. 

We chose the Riccati-Bernoulli equation because it illustrates the process so perfectly 
and because a direct estimate of the accuracy of the partial sum of the series solution is 
troublesome. 

We did not investigate the accuracy for the negative half of the interval, i.e., for —0.2 ^ 
X ^ 0, which we leave as an exercise for the reader. The only change occurs in the estimate 
of y"{x) since then x can be equal to negative numbers. 

Finally we observe that we did not exploit the sign of the error. In fact, our estimate 
(I5.3ip allows us to say that the error we commit is between an error in defect smaller than 2 
units in the eleventh decimal place and an error in excess smaller than 8 units in the twelfth 
decimal place. Therefore we can centralize the error by adding the term 

281 - 665.9 10 _ -1283 
2 ■ 10! ^ ~ 24192000^^ 

to our polynomial (15.320 and obtain an approximating polynomial whose maximum error is 
between ±1.4 units in the 11th decimal place. 
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